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Fractal Structure of Equipotential Curves on a Continuum Percolation Model 
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We numerically investigate the electric potential distribution over a two-dimensional continuum 
percolation model between the electrodes. The model consists of overlapped conductive particles on 
the background with an infinitesimal conductivity. Using the finite difference method, we solve the 
generalized Laplace equation and show that in the potential distribution, there appear the quasi- 
equipotential clusters which approximately and locally have the same values like steps and stairs. 
Since the quasi-equipotential clusters has the fractal structure, we compute the fractal dimension 
of equipotential curves and its dependence on the volume fraction over [0, 1]. The fractal dimension 
in [1.00, 1.257] has a peak at the percolation threshold p c . 
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I. INTRODUCTION 

In the series of articles [2, 0j , we have studied the elec- 
tric conductivity of a percolation model using the finite 
difference method (FDM) [3| ■ The purpose of these stud- 
ies is to reveal the electric properties over the continuum 
percolation model (CPM). As we explored three dimen- 
sional shape effects on the conductivity on CPMs in the 
previous articles [J [2| , in this article we investigate the 
fractal properties of the the solutions of the two dimen- 
sional generalized Laplace equation, 



V • crV0 = 0, 



(1) 



for the conductivity distribution a(x, y) of CPM. Re- 
cently materials consisting of the conductive particles in a 
base material with extremely high resistance are studied 
[J-Q- We handle the binary local conductivities a(x, y) of 
the binary materials in CPM, i.e., conductive particles on 
the background with an infinitesimal conductivity. The 
infinitesimal conductivity i.e., the extremely high resis- 
tance, avoids indeterminacy of the solutions of ([T]) all 
over there, because approaching to zero differs from the 
zero itself. 

Two dimensional CPMs of the overlapped circles were 
studied well in Refs0 and § and references therein. In 
this article, we also deal only with CPM whose particles 
are overlapped circles with the same radius. 

Since the solution of the generalized Laplace equation 
([I} provides the electric potential distribution, we show 
that there exist quasi equipotential clusters, in which the 
potential behaves like steps and stairs approximately. 
There are many pieces whose potentials are approxi- 
mately flat because the electric connections are very 
small and thus the resistance is so high at the bound- 
ary whereas the inner resistance of the particle or the 
composite particles are much small. 
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The existences of the quasi-equipotential clusters 
means the anomalous behaviors of the equipotential 
curves. Though the relation between the percolation 
cluster and fractal structure has been studied Q , in this 
article, we provide the results how the fractal dimension 
of the equipotential curves behaves depending on the vol- 
ume fraction p G [0, 1]. Under the percolation threshold 
Pc our computation can read the computation of dielec- 
tric behavior on a random configuration of metal parti- 
cles in the dielectric matter, which was reported in Refs.[9| 
and [l0| if we interpret the conductivity as the dielectric 
constant. 

Contents in this article are as follows. Sec. II shows 
our computational method; Sec. II. A is on the geomet- 
rical setting in CPMs and Sec. II. B gives the compu- 
tational method of the conductivity over CPMs using 
FDM, which is basically the same as that in the previous 
article |l|. Sec. Ill shows our computational results of 
the solutions of (JXJ) in CPMs and discuss our results from 
physical viewpoints. 



II. COMPUTATIONAL METHOD 

In this section, we explain our computational method 
of CPM. Though we deal only with two dimensional case, 
its essential is the same as the three dimensional case 
using FDM [3f whose detail is in the previous article [l|. 



A. Geometrical setting 

We set particles parameterized by their positions (x, y) 
into a box-region B :— [0, xq] x [0, yo] at random and get 
a configuration 9\ n as one of CPMs. In this article, we 
set xq = yo = L = 100. The particle corresponds to 
a stuffed circle with the same radius p = 1, B XitVi := 
{(x,y) e B | \{x,y) - (Xi,yi)\ < p}. The configuration 
9\ n is given by 9\ n := Ui=i Bxt,vn where each center 
(xi,yi) is given at uniform random in B. We allow their 
overlapping. 
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FIG. 1: The configurations of the seed i s = 1 with the volume 
fractions 0.5, 0.55, 0.6, 0.65, 0.7, 0.8 for (a), (b), (c), (d), (e), 
(f) respectively. 



By monitoring the total volume fraction which is a 
function of 9\ n and is denoted by Vol(fH n ), we continue 
to put the particles as long as Vo\(9\ n ) < p for the 
given volume fraction p. We find the step n(p) such that 
Vol(JH n ( p )_i) < p and Vol(*H„( p )) > p. Since the dif- 
ference between Vol(9^ n ( p )_ 1 ) and Vol(9i n (p)) is at most 
3 x 10~ 4 for the employed parameters. Thus we regard 
Vol(9l„(p)) as each p hereafter under this accuracy. 

Since we use the pseudo-randomness to simulate the 
random configuration 9\ n ( p ) for given p, the configuration 
^n(p) depends upon the seed i s of the pseudo-random 
which we choose and thus let it be denoted by d\ P: i s . For 
the same seed i s of the pseudo-random a configuration 
9l Pj i s of a volume fraction p naturally contains a configu- 
ration 9tp',i s of p 1 < p due to our algorithm. Hence the el- 
ements in the set of the configurations {9lp.i s P S [0, 1]} 
with the same seed i a are relevant. Fig. 1 illustrates the 
configurations of the seed i s = 1 for several p. 



FIG. 2: Conductivity curves; (a) for the linear total conduc- 
tivity atotai and (b) for its logarithm scale. 

in Ref[l|, we numerically solve the generalized Laplace 
equations (TTJ) for the conductivity distribution <j{x, y) for 
each Dip,i s ■ 

Then we obtain the total conductivity Ctotai of the sys- 
tem after we integrate the current <rV0 over the line par- 
allel to the x-axis. Since Ctotai is determined for each 
9^p,i B > Ctotai is a function of the volume fraction p and the 
seed i s of the pseudo-random. We denote it by cr to tai(p)- 

The dependence of the total conductivities on the vol- 
ume fraction as a conductivity curve is given by 



CTtotal(p) = 



(1.0- 

o 



Pc 



for p e [p c , 1], 
otherwise, 



(2) 



where the p c is the threshold and t is the critical expo- 
nent, or merely exponent. Since the threshold p c and the 
exponent t are determined by {£Kp,i s | p £ [0,1]}, they 
are also functions of the volume fraction p and the seed 
i s of the pseudo-random. 



B. Computation of the potentials in CPM 

To apply FDM to the generalized Laplace equation in 
CPM, we use a 2000 x 2000 lattice to represent the square- 
region B. The radius of the particle p = 1 corresponds 
to 25 meshes. 

Further we set the binary local conductivities cr(x, y) 
which consists of the conductive particles with the con- 
ductivity density <7 ma t = 1, and the background with 
infinitesimal conductivity crj n f = 10~ 4 . 

In order to compute the potential distribution ^(cc, y), 
we set 0=1 and ip — on the upper and the lower 
segments respectively as the boundary condition corre- 
sponding to the electrodes. As the side boundary condi- 
tion, we use the natural boundary for each y-boundary 
so that at the side boundaries, current normal to each 
boundary vanish. 

Following our algorithm of FDM as mentioned in detail 



III. THE RESULTS AND DISCUSSIONS 

A. Conductivity curve 

The ordinary computations of the conductivity over 
the percolation system are performed using the Y-A 
algorithm [ll| by assuming that the the conductivity of 
the background insulator vanishes [12|, in other words, 
([T]) is not defined over there in the ordinary approaches. 

In our computations using FDM, we handle the binary 
conductivities <j{x.y) due to the demands of the recent 
material science j^-l6|- 

Even though we deal with the binary conductivities 
with finitely different resistances, the conductivity curve 
is naturally obtained as in Fig. 2 and obeys the (UJ which 
exhibits the properties of the percolation system, because 
the ratio between the binary conductivities <7 ma t and <7; n f 
is sufficiently large. Fig. 2 (a) exhibits the linear scale be- 
havior of the total conductivity whereas Fig. 2 (b) shows 
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FIG. 3: Potential distributions of the seed is = 1 with the vol- 
ume fractions 0.5, 0.55, 0.6, 0.65, 0.7, 0.8 for (a), (b), (c), (d), 
(e), (f) respectively. The black corresponds to 4> = whereas 
the white does to cf> = 1. The graduated gray interpolates 
them. The quasi-potential clusters exist there. 



FIG. 4: Equipotential curves of the potential distributions of 
the seed i 3 = 1 with the volume fractions 0.5, 0.55, 0.6, 0.65, 
0.7, 0.8 for (a), (b), (c), (d), (e), (f) respectively. The curves 
correspond to the values <j) = 0.1, 0.2, . . ., 0.8, and 0.9. 



the logarithm properties. Even though the Fig. 2 (b) il- 
lustrate the properties of the binary materials, the linear 
scale behavior of the total conductivity is described well 
by (|2|). We computed twenty cases with different seeds 
i s of the pseudo-randomness and determine the thresh- 
old p c and the exponent t using the least mean square 
method for the points p > 0.6. For the simpicity, in this 
article, we display only the first five cases as in Fig. 2 
because they represent the others well. 

The average of the percolation threshold is 0.66 ± 0.01 
which slightly differs from 0.674 ± 0.04, 0.688 ± 0.005^ 
0.676 ±0.002, 0677 and 0.691 which are reported in RefQ 
and the exponent is given as 1.33, which recovers the uni- 
versal value 4/3 of two-dimensional percolation [7| with 
this accuracy. Since the effect of the finite size of the an- 
alyzed region makes the threshold smaller, the difference 
comes from the finite size effect |l|, [2j but is not essential. 



B. Potential distribution 

By numerically solving ([1]), we obtain the potential 
distribution <f)(x, y) for each volume fraction p 6 [0, 1] 
and the seed of the pseudo-randomness as we display the 
results in Fig. 3. Fig. 4 shows its equipotential curves for 
= 0. 1,0.2,..., 0.9. 

Fig. 5 is parts of Fig. 4 with the configurations of par- 
ticles at p — 0.5 < p c and p = 0.8 > p c . The contours 
basically run to avoid the conductive particles due to the 
difference of the conductivity. Especially in Fig. 5 (a) for 
the case p = 0.5, the avoidance explicitly appears. 

On the other hand, at p — 1, the potential is simply 
described by 4>(x, y) = y/yo and the contours penetrate 
into the conductive particles. Thus even for p = 0.8, 
the contours penetrate into the conductive particles on 
the conductive paths i.e., the percolation clusters which 
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FIG. 5: The equipotential curves with 8(f) — 6.01, and the 
configurations of particles; (a) for p — 0.5 and (b) for p — 0.8. 



are connected with both anode and cathode. However 
around the isolated percolation cluster which is not con- 
nected with electrodes, the contours go through the gaps 
among the conductive particles. Fig. 5 (b) contains the 
similar states to those situations, in which the contours 
partially penetrate into the conductive particles and par- 
tially avoid the conductive particles. However due to the 
above considerations, these behaviors in Fig. 5 (b) can 
be naturally interpreted. 

Further Fig. 5 shows that the local equipotential curves 
depends upon the local configurations of these particles. 
For given resolution and local region, the sets of particles 
surrounding the region locally play the role of the elec- 
trodes but the direction of the locally averaged gradient, 
i.e., the local electric field, strongly depends upon the 
resolutions and upon the region which we average. Since 
XY-model which is the 0(2) valued spin system also have 
the fractal structure [13J , the randomness brings the self- 
similar structure into our system. 
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FIG. 6: The fractal dimension v.s. volume fraction p for five 
seeds. 



As Mandelbrot showed that the fractal picture is re- 
lated to geographical geometry in Ref. [lj, Fig. 3 looks 
like pictures of the piled up mountains in distance using 
the monochrome water-painting. There are many pieces 
whose potentials are approximately flat, like leaves of the 
piled leaves. These properties are interpreted that there 
exist quasi-equipotential clusters whose individual poten- 
tial is approximately same and behaves like steps and 
stairs. In other words, the electric connections among 
the pieces are very small and thus the resistance is so 
high at the connections whereas the inner resistance of 
the particle or the composite particles is much small. 



C. Fractal dimension of equipotential curves 

Let us consider the effect of the quasi-equipotential 
clusters on the equipotential curves in Figs. 4 and 5. 
Under the threshold p c , this phenomenon is essentially 



the same as the one of dielectric properties related to the 
electric breakdown, which were studied by Gryure and 
Beale in Ref. [9| and recently by Stoyanov, Mc Carthy, 
Kollosche and Kofod in Ref. JJLQfc though the authors han- 
dled the hard core model in which the overlapping is 
prohibited instead of our soft core CPMs. They stud- 
ied the dielectric materials as background with infinite 
dielectric particles (metal particles), which are governed 
by the same generalized Laplace equation ([1]), though it 
is for the region of p € [0,p c ]. Fig. 5 (a) corresponds 
to the Fig. I. in Ref. [9| if the conductivity reads the 
dielectric constant. 

In our computational condition on the conductivity, 
we can determine the potential distributions over B for 
every volume fraction p <E [0, 1] due to the infinitesimal 
conductivity <7i n f. 

As objects in the percolation theory exhibit the fractal 
structures, related to the percolation clusters [8{, related 
to the breakdown [L5l - fl8j and so on, the fractal dimension 
characterizes these properties. We also computed the 
fractal dimensions but especially ones of the equipotential 
curves. We used the box-counting method [14|. Since 
the curves in Fig. 4 of the seed i s — 1 have the fractal 
structure, their fractal dimensions are not trivial. 

The dependences of the fractal dimensions of the 
equipotential curves on the volume fraction p £ [0, 1] for 
five different seeds i s are presented in Fig. 6. Mandel- 
brot showed that the geographical curves like the coast- 
lines sometimes behave fractal [14| and have the fractal 
dimensions 1.0^1.25. Our fractal dimension is in [1.00, 
1.257] and has a peak at the threshold p c . In other words, 
around p c , the equipotential curves behave complicatedly 
and have the self-similar structures. On the other hand, 
for the p = and p = 1 cases, the dimension must be 
equal to 1. Thus the shape of the graph in Fig. 6 shows 
the curve with a peak at the percolation threshold p = p c . 
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